Pseudo-Hermitian version of the charged 
harmonic oscillator and its "forgotten" exact 

solutions 



' Miloslav Znojil 

o i 

Nuclear Physics Institute of Academy of Sciences of the Czech Republic, 



O 



£^ ' 250 68 Rez, Czech Republic 

e-mail: znojil@ujf.cas.cz 

m 



> 
in 
oo 
o , 

Abstract 

An unusual type of the exact solvability is reported. It is exemplified by the Coulomb 
plus harmonic oscillator in D dimensions after a complexification of its Hamiltonian 
which keeps the energies real. Infinitely many bound states are found in closed form 
which generalizes the popular harmonic-oscillator states at zero charge and even 
parity. Apparently, the model is halfway between exact and quasi-exact. 
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1 Introduction 



Schrodinger equation for the shifted and charged harmonic oscillator in D dimensions 
reads 



d 2 £{£ + !) f n 
T1+ 9 + - + 2gr + r 
ar r z r 



ip(r) = E ip(r), re (0,oo). (1) 



Thirty years ago, Andre Hautot ||T| noticed that at certain non-vanishing couplings 
/ and/or g it may posses elementary solutions for the equidistant set of the energy 
levels 

E = E f; utot) = 2n + 2 + C-g 2 , (2) 

n = 0, 1, . . . , £ = £{£) = 2£ + 1 = D - 2, D, D + 2, . . . 

which do not depend on the charge /. The charge itself is not arbitrary (for this 
reason, the models of this type are called quasi-exactly solvable, QES). One has to 
evaluate the admissible values of / = f n as real roots of certain polynomials of the 
(n — 1)— st degree (see below). 

The undeniable mathematical as well as physical appeal of QES solutions has 
been revealed by many independent authors whose work was summarized by Alex 
Ushveridze 0. Very recently, the QES models helped to clarify some counterintuitive 
formal features of the so called VT symmetric quantum mechanics of Bender et al 
who replaced the usual Hermiticity H = W by the mere commutativity of the 
Hamiltonian with the product of parity V and time reversal T. In the early stages of 
development, the studies of this formalism were strongly motivated by its relevance in 
field theory |4j]. In such a setting, it was very impressive when Bender and Boettcher 
demonstrated that in contrast to the current Hermitian case, quartic polynomial 
oscillators belong to the QES class after their appropriate VT symmetrization (cf. 
also ref. |J for more details). 

The charged harmonic oscillator ([!]) does not possess similar appeal in field theory 
but it was still amusing to reveal in paper that its Schrodinger equation does not 
lose its partial elementary solvability even after the weakening of the Hermiticity 
to the mere VT symmetry of its Hamiltonian. In accord with the Buslaev's and 
Grecchi's recipe M we used the shifted coordinates x 6 (—00,00) on complex line 
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r(x) = x — i e at a constant distance e > from the real axis 
d 2 £(£+1) . R 



, 9 H tt^ - + i -rr + 2ibr(x) + r (x; 

ax" 1 r^(x) r(a;j 



i> n (x) = E n i/j n (x). (3) 



Unfortunately, we only analyzed the £ = solutions in the quasi-odd regime (cf. a 
more detailed explanation below). 

Since the publication of paper a significant progress has been achieved in the 
interpretation of the non-Hermitian equations. Several authors || emphasized that 
within the domain of quantum mechanics, the VT symmetry of eq. (|3|) should be 
replaced by the (formally equivalent but mathematically more natural) requirement 
of the pseudo-Hermiticity of the Hamiltonian, 

H = rj^H^T), r] = rj ] . (4) 

One of the oldest illustrations of the efficiency of the use of the pseudo-Hermitian 
Hamiltonians (|j) with an auxiliary indefinite metric r\ appears in the relativistic 
quantum mechanics where it very naturally originates from the Feshbach's and Vil- 
lars' Hamiltonian formulation of the Klein Gordon equation for the particles with 
zero spin |H|. The same choice of a suitable invertible rj helps to clarify some con- 
temporary problems in the cosmological models based on the equations of Bryce de 
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In the light of eq. (|J), physical interpretation of the non-Hermitian bound states is 
more transparent and does not depend too much on the specification of the operator 
rj itself. This operator only plays the role of a certain auxiliary transformation of the 
dual Hilbert space. For a detailed explanation of this idea we recommend the older 
review |12| where the physical meaning of the nontrivial "metric" rj was illustrated 
by its emergence in the many-fermion models where rj ^ I characterizes the so called 
Dyson's mapping of the "physical" (and Hermitian) fermionic Hamiltonians onto 



their "more easily solvable" r]— Hermitian bosonic equivalents [|T3 

One must keep in mind the non-uniqueness of the metric 77 which belongs to the 



given H. According to A. Mostafazadeh |L4| one can replace the initial indefinite 



(and, in particular, VT— symmetric) metric rji by an alternative r] 2 which is posi- 
tive definite. In the other words, the puzzling quasi-unitary evolution generated by 
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the indefinite 771 ||T5| may be declared an artifact of our constructions. Vice versa, 
all the phenomenological considerations should necessarily be related to the positive 
definite version r]2 > of the metric (one can show that is exists for all the diago- 



nalizable Hamiltonians H |14}l). Then, the time evolution remains compatible with 



the probabilistic interpretation of the norm in the Hilbert space of states. 

In the light of the possible peaceful coexistence of rji = V with some 772 > in 
eq. our attention has been re-directed to the solutions of eq. which diverge in 
the simple-minded Hermitian limit e — > and which were omitted from the scope of 
our preceding study JF]. We are going to correct the omission now. Our expository 
section ^] will outline an improvement of the method which generates the old QES 
solutions of ref. [[7]. Section [3] will then modify the basic ansatz which opens the 
way towards the new (or, in the wording of our title, "forgotten") QES solutions of 
eq. @. Some of their properties and possible applications will be finally discussed 
in section f|. 



2 The standard quasi-exact solutions 

As long as the differential eq. (|3|) is of the second order, its general solution is a 
superposition of some two linearly independent solutions. This independence may 
be deduced from their available leading-order form near the origin, 

^( r ) = c.^-\r) + c + ^+\r), ^±\r) = O ( r Va=F(*+i/2)) . (5) 

In the spirit of ref. |7j one usually searches for the polynomial solutions compatible 
with the correct physical boundary conditions in the origin (i.e., c + = |16|]) as well 
as with their asymptotic normalizability. In the spirit of the general QES philosophy 
one may meet both these requirements by employing the special elementary ansatz 
of the very common harmonic-oscillator-like form 

N 

^(x) =^-\x) =e~ r2/2 ~ ibr hi n\ ir T +e+ \ r = r(x)=x-ie. (6) 

n=0 

The construction of the solutions of this type degenerates to the insertion of eq. (0) 
into the differential eq. (§) which gives the homogeneous set of N + 2 linear algebraic 
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equations for the N + 1 coefficients hi \ We may drop the superscript and turn to 
the explicit non-square matrix form of the latter equations, 



1 B C 
A, B l 



\ 



A 



N-l 



Bn-i 
A n 



V 



Cn-i 
B N 

A N+ i ) 



u \ 

h 2 



\ h N j 



(7) 



with elements 

A n = = b 2 + 2n + £ - E, B n = flH 
C n = CH = (n + 1) (n + 1 + C), L = 2l- 



-- -(2n + 1 + C)b - F, 
1, n = 0, 1,... . 



This is a finite-dimensional and over-determined linear algebraic re-incarnation of 
the original differential equation (|3|). Its matrix structure enables us to define the 
wave function (i.e., all its energy- dependent coefficients) as determinants, 

I Bi n \ 



h 



N 



bj-1 



{-A j )(-A j+1 )...(-A 



det 



N, 



A j+1 B j+ i C j+1 



V 



An-i Bn-i Cn-i 

An B n J 

with j — N, N — 1, . . . , 1 and under any choice of the normalization hN 7^ 0. 

The latter normalization convention converts the last row An+iIin = of eq. (|7|) 
into the constraint 

E = E^ = 2N + 2 + C + b 2 . (9) 

In the other words, the condition of the mutual compatibility of the original over- 
determined linear system (0) fixes the energy which coincides with the old Hautot's 
formula @. At any N = 0,1, . . . the energy is an increasing function of the angular 
momentum i or £ and of the size of the shift b. The QES construction is complete 
and 

• simplifies the prescription of ref. (where the special cases of eqs. (0) and 
(PD contained the less compact matrix with four diagonals), 
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• leads to the polynomial wave functions (with the closed form (p|) of the coeffi- 
cients), 

• preserves the Hautot's explicit form (|14D of the energies, 

• reduces the differential Schrodinger equation to its n— dimensional square- 
matrix form. 

The first observation (simplification) is a marginal technical merit due to our tran- 
sition to a better ansatz. In contrast, the last feature of the QES solutions remains 
highly unpleasant as it forces us to guarantee that the related secular determinant 
vanishes, 



i(-) 



det 



a 



(-) 

D 

(-) 



A[~> Bf C{ 



y(-) 



A 



(+) 

N-l 



B 



(-) 

N-l 



C 



(-) 



V 



.4 



(-) 

N 



B 



N-l 
(-) 



N 



. 



(10) 



/ 



Such a constraint determines the set of the N + 1 admissible couplings F = F% ■* (N) 
Jl| and its purely numerical nature is an example of the most serious practical short- 
coming of the majority of the QES models 0. We are now going to describe a 
remarkable exception from this discouraging rule. 



3 Nonstandard, quasi-even QES states 
3.1 The concept of quasi-parity 

Above we emphasized that at a fixed, non-vanishing shift e > the ambiguity of the 
metric rj opens the possibility of using the "simpler" conjugation ([|) with r\ — r\\ — V 
during the explicit constructions of the solutions while switching to their "physical" 
re-interpretation based on an alternative scalar product with the positive definite 
metric r] 2 > 0. In such an approach both components c_ 7^ and c + 7^ of wave 
functions in eq. (||) may be equally useful. 

Once we relax the redundant boundary conditions in the origin we get more 
solutions of course. One of the most transparent illustrations of the emergence of 
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the additional "quasi-even" c + 7^ solutions was described in ref. [|17j where eq. (|3]) 
has been solved at the vanishing b = F n = 0. The superscripts in the resulting 
states ?/4 + )(r) and ip^~\r) have been interpreted as the so called quasi-parity. The 
introduction of this concept was motivated by the observation that the quasi-parity 
degenerated to the current parity at b = F n = £ = (cf. eq. (0)). 

Its independent additional support appeared in ref. [0 revealing its connection 
with the VT parity and with the pseudo-norm using rj = rji = V. It even plays its 



role in the supersymmetric quantum mechanics (cf. ref. |i| for more details) but is 
missing from our ansatz (||) inherited from ref. 0. As we understand it now, ansatz 
([]) is unnecessarily restrictive as it represents merely quasi-odd solutions. From its 
generalization 

N 

Mr) = e- r2 / 2 ~ tbr J2 h n {ir) n ~ l (11) 

n=0 

(where we dropped all the superscripts for the time being) one can always return to 
the old quasi-odd option via the additional superscripted constraint 

4~ } = h[- ] = ... = h { c \ = . (12) 

Whenever necessary, the superscripted "quasi-even" QES solutions may be char- 
acterized by the alternative criterion 

+ + + >0. (13) 

We are now close to our key claim that the structure of the quasi-even QES solutions 
is exceptionally simple. 



3.2 QES states having the even quasi-parity 

The source of the latter claim lies in the improved ansatz ([□]) which leads to the 
same equation (|7|) with the new matrix elements 

A n = = b 2 + 2n - C - E, B n = 5« = -(2n + 1 — C)b — F, 
C n = CW = (n + 1) (n+ 1 - £), C = 2£+l, n = 0, 1, . . . . 

The energy formula @ is only marginally modified, 

E = E {+) = 2N + 2- C + b 2 . (14) 
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Still, we immediately notice the much more important difference connected with the 
presence of the vanishing matrix element = in the upper diagonal of our new 
form of the QES secular equation. This means that the ^—superscripted secular 
determinant may be re-written as the product of a "small" , /^—dimensional 



S {s) = det 



4 +) b[ +) c\ 



(+) 



A 



(+) 

C-2 



\ 

with another determinant of a "large" 
factor 



B 
A 



(+) 

C-2 

(+) 
C-l 



C 



(+) 

C-2 



B {+) i 



( Bi +) 



S {L) = det 



.4 



(+) 
c+i 



C 
B 



(+) 
c 

(+) 

C+l 



(N — £)— dimensional matrix. The latter 
\ 



r (+) 



A (+) 



B 



(+) 



N-l 



N-l 



V 



A 



(+) 

N 



B 



(+) 

N 



J 



precisely coincides with the left-hand side expression in eq. flT0| ) which guarantees, 
in its turn, the existence of the QES solutions with the property (0) . The use of the 
condition = would return us back to the old quasi-odd ansatzs of section 0. In 
what follows we shall ignore these solutions as standard and assume that ^ o. 



3.3 Facilitated QES constructions 

After one concentrates attention solely to the quasi-even states, the QES construction 
degenerates, basically, to the secular equation <S^ = 0. One verifies easily that 
the acceptance of this condition is consistent with the quasi-parity ([131) . The wave 
function coefficients themselves remain formally determined by the ^ —superscripted 
version of the determinants (^) whose dimension grows with N. For this reason we 
recommend a re-interpretation of these coefficients as quantities evaluated by the 
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recurrences initiated at an initial hffi ^ and defining /i)v-i, h^i 2 - ■ ■ ■ s ^ e P by step 



(+) !,(+) 



-2JV /3x-F 4-2£ 



\ 



(+) 



h N _i 

V 4 +) / 



o . 



(15) 



_4 (5 N ^-F N 2 -NC 
y -2 /3jv-F y 

We abbreviated here /3 n = — (2n + 1 — C) b and note that in the last step of these 
recurrences, the secular equation = replaces the redundant definition of a 
ghost coefficient h_i = 0. Finally, the explicit form of our secular equation 



( Po-F l-C 
-2N Pi-F 



det 



4- 2C 



V 



-2{N + 3-C) Pc-2 -F l-C 

-2(N + 2-£) Pc-i-F ) 



(16) 



specifies the family of the admissible charges Fn = F^ (N) in a way which remains 
virtually purely non-numerical for the first few dimensions L. 

3.3.1 C = 2 

One of the key reasons why both the formal appeal and practical importance of the 
quasi-even spectrum (|I4T) remained unnoticed in ref. was purely psychological. 
Indeed, at the simplest choice of C = 1 in eq. flTEl ) (which may mean both the s— wave 
in three dimensions and an even state at D — 0) one does not obtain anything new. 
Equation (16) provides the single root Fq (N) = and one just reveals the well 
known fact that at the vanishing eigencharge our model degenerates to the linear 
harmonic oscillator defined on the whole line. 

Let us move, therefore, to the first nontrivial choice of £ = 2 corresponding to 
the p— wave in two dimensions or to the s— wave in four dimensions. This gives the 
following two series of the fully non-numerical eigencharges, 



F [h2] (N) = ±J(b 2 + 2N). 



(17) 
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The wave functions retain the even quasi-parity in a way compatible with eq. 



h$l 2] (N) = ~ [b + F M (N)\ h^> >2] (N). 

Both the eigencharges grow with the increasing size of the shift b and with the number 
iV of nodes (i.e., with the energy E N ). In the coupling-energy plane our QES states 
may be visualized as families located along certain curved lines which, in a way, lie 
somewhere in between the harmonic-oscillator F = straight line and the numerous 
Hautot's sets of roots Fk(N) each of which is defined at a fixed energy or integer N. 
Sometimes, the similar families of the bound states with the "energy = constant" 
property are being called Sturmians. In this sense one could speak here about a 
certain further generalization of the latter concept. 

3.3.2 C = 3 

A mild formal shortcoming of the present Coulomb + harmonic oscillators lies in 
a quick growth of complexity of eq. ( Jl6| ) for the larger C At any C > 3 one 
should not be tempted to generate the formulae verifying, say, the smooth iV— and 
b— dependence of the eigencharges. In practice, the other approaches may definitely 
prove preferable. 

Even at the very first L = 3 the comparatively compact form of our eq. fll6|), 

( 2b-F -2 \ 

det -2N -F -2 =0, (18) 

v -2N + 2 -26 - F j 

should not inspire a search for the triplet of charges {F , F\, F 2 } via the closed (i.e., 
Cardano) formulae (the reader is recommended to try to generate them using the 
computer symbolic manipulations in order to see that they really are enormously 
clumsy). A significantly better strategy consists in an elimination of the (unique) 
value of N from the above secular determinant (|j|) giving 

N = N(F, b) = -— (AFb 2 + 8b - F 3 - 4F) . 

After we fix any left-hand-side integer iV we may pick up the two eigenshifts frri^i as 
elementary functions of the indeterminate variable F (we skip the details which are 
trivial). 



1 



,(+) 
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3.3.3 £ > 3 

The main advantage of the semi-implicit techniques of the solution of eq. (|16|) is that 
they may work at a few larger integers £ > 3. For illustration, the choice of £ = 5 
(which corresponds to the d— wave in three dimensions) leads already to the purely 
numerical determination of eigencharges F. In an alternative approach we eliminate 

-7686 - 256F6 2 + 768F + 40F 3 ± 24^(1024^ + 1926F 3 + 512F 2 + F 6 ) 

and recommend the graphical determination of the eigenvalues F = F(b) and/or 
b = b(F) afterwards. 

The most practical possibility consists in a direct selection of a suitable shift b 
followed by the subsequent diagonalization of the purely numerical matrix. In an 
illustration using £ = 3 and b = 5 one gets the three eigencharges 

F = {10.757, -10.400, -0.35755} 

at the smallest possible N = 2. In an opposite-extreme using N = 1000 the computed 
values 

F = {89.98, -89.975, -0.0049407} 
already lie very close to their large— N estimate obtainable in closed form, 

F « {V8N, -V8N, -b/N} w {89.44, -89.44, -0.005}. 

This simplifies the verification of the reality of the eigencharges and confirms the 
smoothness of their iV— dependence. Such a type of calculation is very quick and 
gives results sampled at £ = 4 and b = 5 in Table 1. 

4 Discussion 

The main merit of our key eq. ( |16|) (which defines the eigencharges) is that its 
dimension is independent of the quantum number N. Equation dl7f) is the best 
illustration of the related new form of the solvability which we intended to describe 
here. Still, the principle of the whole construction is more general and one might 
try to apply the similar recipe to the quartic oscillator of refs. [El, to the sextic 



512F 
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oscillators studied by many authors [[HJ], to the decadic oscillator of ref. |21[ and to 
the numerous existing modifications of these most popular or "canonical" models. 



4.1 Non-orthogonal QES states as a basis? 

In our particular example, the strength of the Coulombic interaction appears to be 
an energy- or N— dependent quantity. One deals with an L > 2 generalization of 
the common £ = 1 harmonic oscillator. Its most important features are an apparent 
completeness "in a relevant subspace" (a guess inspired by their infinite number) and 
a compact form (reflecting the N— independent evaluation of the eligible charges, 
each of which is selected as a function of the main quantum number N = 0, 1, . . .). 
Both these features make our infinite set(s) of the quasi-even states very similar to 
the ordinary harmonic oscillator basis. In the context of concluding remarks, let us 
pay some attention to the possible analogies of the latter type. 

Firstly, due to the non-Hermiticity of the Hamiltonian H(F) = H(0) + F W we 
have to distinguish between the left (= double-ket) and right (= single-ket) QES 
eigenstates, 

[H(0) + F N W] \N) = E N \N), (19) 

((N\ [H(0) + F N W} = E N ((N\. (20) 

The integer N numbers the energies E^ as well as the selected charges F^ = 
F\k]( + '(N) so that the left and right eigenstates exist at the common energies (|14|) 
and charges [say, flTT])]. The wave functions are defined in closed form, as polynomi- 
als in the coordinates (cf. (0)) and in the couplings and quantum numbers (cf. (|8|)). 
This is of paramount importance, making our quasi-even QES solutions extremely 
similar to the even bound states of the exactly solvable chargeless oscillators which 
form one of the most popular complete bases in L 2 (0, oo). 

In a tentative parallel one could search for the appropriately weakened biorthog- 
onality relations. This is a real mathematical challenge since our QES solutions are 
only defined at the exceptional and I— and N— dependent charges. Still, one can 
easily verify the manifest non-orthogonality of the pairs of many randomly selected 
QES states. Of course, the closed form as well as the infinite number of these states 
inspires their use in a perturbative or variational context. 
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For the similar purposes one has to truncate their set to a finite subset (A/" < oo), 
assuming that the matrix of their overlaps 

Qm.n = ((m\n), m,n = 0, 1, . . .J\f 

is invertible. We may then proceed, say, in the variational spirit and reduce our 
Hilbert space and its dual to the finite-dimensional subspaces spanned by our subset 
of the selected QES eigenvectors. The approximate identity operator becomes defined 
by the usual series 

N 

I = \ m ) R m,n{{n\ 

m,n=Q 

where R = Q~ l is, in general, fully non-diagonal. Also the Hamiltonian H(F) itself 
becomes approximated by a non-diagonal matrix. At almost all F, the search for 
the energies E = E(F) becomes, therefore, a numerical task. 

4.2 Speculations about applicability 

In spite of the unpleasant character of the latter conclusion, one should still feel the 
difference between a fully general matrix diagonalization and our "next-to-solvable" 
Coulomb + harmonic problem considered at any charge F, 

[H(0)+FW] \V) = E(F)\V) . (21) 

In particular, our Schrodinger equation may be (e.g., perturbatively) connected to 
its special cases with QES character. In an attempted step towards making such a 
connection explicit, let us imagine that equations ( |T9D and ( |20| ) share the energy and 
charge (though not the eigevectors) at every given N. This means that, respectively, 
we have the relations 

((N\ [H(0) + F M W] \M) = E M ((N\M), (22) 

«iV| [H(0) + F N W\ \M) = E N ((N\M) (23) 
the subtraction of which gives the strongest constraint 

(F M - F N ) ((N\ W \M) = (E M - E N ) Q N , M . (24) 
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This relation is an immediate generalization of the bi-orthogonality of the states 
which would result from it in the hypothetical case of the subscript-independent 
charges F M = F N . 

Let us now return to the relations (p2|) and ( ^3|) and deduce the matrix form of 
the Coulomb + harmonic Hamiltonian H(F) at any value of the charge, 

((N\ H(F) \M) = (F - F M ) ((N\ W \M) + E M Qn,m • (25) 

This formula will help us to avoid the tedious and most time consuming part of 
the diagonalization of H(F), viz, the evaluation of the input matrix elements which 
precedes the solution of the problem fl2"T|) studied in its explicit matrix form which 
defines, in principle, all the necessary components p^ = ({N\*f>) of the wave functions 

N 

|$) = ^ \m) R m ,nPn- 
m,n=0 

At this moment we have to re-emphasize that we do not intend to perform any 
numerical calculations. We rather wish to stress the helpful role which can be played 
by the QES states. In such a context, the number of the necessary input matrix 



elements encounters the most drastic reduction after the insertion of eq. (j25 ) in 



eq. (H, 

£ ( F ~ F n) W \K) R k ,j Pj ={E- E n ) p N . (26) 

K,J 

Next, the necessary input information is further reduced by the generalized biorthog- 
onality relations (|24l) which express all the off-diagonal elements of the (in our ex- 
ample, Coulombic) operator W in terms of the known overlap matrix Q. In this way 
the numerical or perturbative diagonalization of the matrix Schrodinger eq. Q2"B]), 



(f-f n )y: 
j 



En — E K 



pj = {E- E N ) p N , (27) 



T/v Rn,j + X/ zt - Qn,k Rk,j 

will require the independent input evaluation of the mere diagonal matrix elements 
T/v = ((N\W \ N) assuming of course that we always have Fm ^ F^ for M ^ N. We 
may summarize that the main merit of the use of the QES states lies in the compact 
and easily generated matrix form of our Coulomb + harmonic matrix Schrodinger 
equation at any non-QES charge F. 
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4.3 Outlook 

We have seen that our system (|3|) may quite efficiently be treated by the non- 
numerical as well as almost purely numerical means. In the former sense it has 
been shown to lie somewhere in between the QES and exactly solvable category (cf. 
Table 2). In comparison with its completely solvable neighbor (i.e., with its own 
harmonic oscillator special case of ref. ]T7[), its nontrivial b ^ and F^O versions 



do not generate all their bound states in the elementary form. Still, in contrast to 
the quasi-odd QES model [0] in Table 2, its present quasi-even partner supports an 
infinite number of bound states in the closed, elementary form. Thus, our example 
(as well as any one of its many analogs) could be assigned a "midway" status in some 
of its applications and interpretations. 

Acknowledgement 

Work supported by GA AS (Czech Republic), grant Nr. A 104 8004. 



14 



Table 1. iV— dependence of eigencharges at C = 4 and 6 = 5. 



N 


Fn 


3 


-15.611 


-5.9279 


4.888 7 


16.651 


30 


-27. 149 


-9.290 9 


8.929 4 


27.511 


100 


-44. 732 


-15 


14. 865 


44. 867 


200 


-61.665 


-20.602 


20.531 


61.736 


300 


-74. 856 


-24. 984 


24. 936 


74. 904 


1000 


-134. 93 


-44. 985 


44. 970 


134. 94 


3000 


-232.82 


-77.610 


77. 605 


232. 83 


30000 


-734. 99 


-245. 00 


245. 00 


734. 99 



Table 2. Solvable pseudo-Hermitian potentials: Tentative classification. 



class 


quasi — exact 


intermediate 


exact 


solutions available 
range of couplings 
illustrative example 


at a finite set of N 
restricted 

E 


at infinitely many iV 
restricted 
here 


at all N 

any 
|17| 
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